Entangling and decohering influence of noisy 
perturbations. Stochastic anti-resonance in the time 
evolution of interacting qubits. 

Demetris P.K. Ghikas*, Athanasios C. Tzemos 

Department of Physics University of Patras, 
Patras 26500, Greece 



Abstract 

We study the entanglement evolution of two coupled qubits in interaction 
with an external environment and in the presence of an external magnetic 
field with a stochastic component. The results show the expected degra- 
dation of entanglement due to the noise. The new effect is that the time 
of disentanglement depends in a non-monotonous way on the strength of 
the noise. We find that it is shortest for an intermediate strength value of 
the latter. This we call "stochastic anti- resonance" . Our results could lead 
to a better undestanding of noisy perturbations and their role for optimal 
designing of quantum devices. 

Keywords: entanglement, Heisenberg chain, noise, stochastic resonance 
1. Introduction 

The source of difficulties for the scalability of quantum information pro- 
cessing devices is the constant interaction with their local environments 
[H, Q, y, 0, S] . Engineering these interactions can improve the performance, 
but residual noisy agents, like non-zero temperature environments, or noisy 
external fields of control, impair the optimality of the behavior. Thus, appart 
from the well established methods of decoupling schemes error correcting 
codes 01 and decoherence free subspaces [sj, addressing directly the role of 
noise is an extra possibility for the engineering methods. The discovery of the 
counterintuitive phenomenon of stochastic resonance, in both classical non- 
linear jij and quantum systems [10], points to the direction of a non trivial 
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role of noise as a (possibly beneficial) control. There are extensive studies of 
entanglement and decoherence in thermal environments and in the presence 
of external magnetic fields and noise [O, 12, 13, 14, 15, 16, 17|. Especially 



the stochastic control and stochastic resonance has been studied in various 

Here we study the dynamical evolution of 



18, 19, 20, 21, 22 



arrangements 

entanglement for a two qubit Heisenberg chain in interaction with an external 
bosonic environment and a noisy external magnetic field. We are interested 
in the evolution of the effects of classical noise on the quantum properties of 
a simple system that can be solved analytically. More specifically, we have 
looked for a possible non-monotonous dependence of some physical quan- 
tity on the strength of an external noise, namely whether we can have a 
stochastic resonance type behavior. The appearance of such a behavior, ap- 
part from its inherent interest, would be useful for a better understanding 
of noisy perturbations in real quantum information devices. What we have 
found is that the noise degrades initial entanglement, something expected, 
but the time taken for this effect is not a monotonous function of the noise 
strength for some initial states. This is interesting since it shows that an 
intermediate strength of noise can be worse than a smaller or bigger. Since 
we observe a worse degradation, we call the effect stochastic antiresonance. 
What in fact we have found is that the Entanglement Sudden Death (ESD) 



time [23|, |2J, |25|, |26|, l27|, |28| does not depend monotonically on the strength 



of the noise. We are currently looking to multiqubit systems with a wider 
range of parameter values and with the use of stochastic differential equa- 
tions for a direct numerical study. In the paper first we introduce the corre- 
sponding master equation where the external gaussian white classical noise 
is incorporated in terms of a double commutator. Then the Heisenberg type 
Hamiltonian is presented with some comments on the parameter values. Our 
results on the concurrence and purity evolution are presented in the form of 
2d and 3d diagrams where the appearance of stochastic antiresonance effect 
is made evident. 



2. Quantum systems driven by external white noise 

In order to be able to have analytic expressions for the dynamics, we use 
a form of the master equation where the noise is incorporated in terms of an 



29j. For a zero temperature bath the master equation 



appropriate term as m 
is 

dt 

where V=Sf + 5|, D[S^]ps = S^psiS^Y - {{S^Y , ps}/2, D[S,]ps = 
S2Ps{S2y — {{.32)^82, Ps}/'2 and 7 is the rate of population relaxation. We 



-t[H, p,] + 7^[5r]p. + lD[S2]ps - M[V, [V, p,]]. 
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assumed that each qubit has the same interaction with the environment. The 
first term describes the unitary evolution of the system, the second and third 
terms the interaction between the system and a thermal environment. The 
last term is the addition of an external classical gaussian white noise [ioj. 

3. Two qubit XY Heisenberg model 

The general form of a N-spin Heisenberg chain (for spin 1/2 particles) 
with nearest-neighbor interaction is 



N 



n=l 



where Sl^ = \cr^{i = x,y,z) the spin-1/2 operators, cr" the corresponding 
Pauli's operators and ^ = 1. Moreover if the chain is periodic, the boundary 
condition is S^^i = Si. The chain is called XYZ model for arbitrary values of 
Jj. If Jj < the chain is called ferromagnetic and if J, > 0, antiferromagnetic. 
Here we study the case where Jz = 0. The Hamiltonian of a two qubit 



Heisenberg XY chain in the presence of an external magnetic field u is |ll 

H = uiS^, + S^) + JiS^S, + Si St) + AiS^St + SiS^), (3) 

where J = (J^ + Jy)/2, A = (J^ - Jy)/2 and S^ = S'^ ± iSy. The first 
term describes the energy levels of spins in the external field, while the other 
two terms are the spin-interaction Hamiltonian and describe the coherence 
production of the two qubits. 



As it has been shown [llj,[15| the third term is essential for the produc- 
tion of a non-vanishing steady-state entanglement even in the presence of 
decoherence. Moreover, one can produce entanglement for any finite tem- 
perature, by adjusting the magnetic field, but only in the anisotropic case 
(A 7^ 0) The XY Heisenberg chain has been studied intensively in the 



last decade, for various purposes, such as NMR quantum computation [30 



and quantum teleportation [16|, because of its simplicity (especially in the 
two-qubit case) and its many interesting features. 

For the sake of simplicity we assumed that 7 is the same for both qubits. 
This assumption requires the use of system parameters, such that the in- 
equalities 



-™'-^<0.1 (4) 
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and 



M<o,i 



(5) 



are satisfied. These equalities ensure that the interaction between the two 
qubits does not alter the energy level separations by more than 10% from 
that of the non-interacting system. ll|] All of our examples satisfy the above 
equalities. 



4. Concurrence and purity evolution for T = 



We use concurrence for the quantification of entanglement [3l[ , [32] , j33j , [34 
For a system density matrix p, the concurrence C is 



C = max(A/Ai — a/A2 — a/As — a/A4, 0), (6) 

where Ai, A2, A3, A4 are the eigenvalues of spin flipped density matrix R (with 
Ai the largest one). The definition is: 

R = p{ay(g)ay)p*{ay®ay) (7) 

C is in the range [0, 1]. C = corresponds to a product state, while C = 1 to 
a maximally entangled state. All the other states inside this range are called 
partially entangled. 

For the purity of quantum states, we use the linear entropy function, 
defined as: 

L = l- Tr{p^) (8) 
The basis states for our system are the usual product states: 



= |e)i$ 


|e)2, 


(9) 


IV'2) = |e)i§ 


^ \9)2, 


(10) 




5 |e)2, 


(11) 


1^4) = \g)i ^ 


5 1^7)2. 


(12) 



A simplifying feature of Hamiltonian ([3j) is that the dynamics of the 
system evolution splits into two independent sets of equations for the two 
submatrices. 
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We note that the second submatrix has not the form of a density matrix. 
In this paper we use the first submatrix for our results, which can describe 
product states, entangled states and mixtures of them. For such a density 
matrix, the concurrence and linear entropy functions are found to be equal 
to: 



C = max{0,Ci,C2}, 



(14) 



where 



C2 



and 



2(|p4i| - VP33P22), 

2(|P32| - VP44P11), 



(15) 
(16) 



L = l-Y,pl{t)-2\p32{t)\'-2\p,,{t)\' 



(17) 



i=l 



correspondingly. 

We present the steady state solution analytically in the appendix. Here 
we report on our results with 2-d and 3-d diagrams corresponding to some 
fixed values of system parameters. For M = we get the results of |ll|. In 
all cases the parameter values are J = A = 0.1,7 = 0.01 and to = l\j We 
have studied the evolution for various values of the noise strength, with initial 
conditions associated with product states. Bell states and mixed states. The 
corresponding representation of the initial density submatrices are as follows. 
Product state \gg)'- 



18) 
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/I 








o\ 


























vo 








0/ 



(19) 



7 <C w so that the weak coupling approximation remains valid. 
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Figure 1: Concurrence as a function of time and noise parameter M for the product state 
= \gg) (left) and Bell state ?/'m = 1^^^^^^ (right) ( J = A = 0.1,7 = 0.01 and cj = 1). 
Noise degrades entanglement monotonously. 



Bell state 



/1/2 l/2\ 





Vl/2 1/2/ 



(20) 



All of the above initial conditions have matrix representation of the form 



/• 





• 00. 



(21) 



(i.e. the evolution begins outside the central subspace {\eg), \ge)}) and do 
not exhibit stochastic antiresonance. Noise degrades entanglement, both the 
initial and final, in a monotonous way. 

For initial conditions with matrix representation of the form 



/O 0\ 

• • 

• • 

\0 0/ 



(22) 
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Figure 2: Concurrence collapse (left) and revival (right) as a function of time and noise 
parameter M for the Bell state V'm — On the left it is clear that there is a noise 

parameter value near 0.5 for which the collapse time is fastest. On the right it is obvious 
that the bigger the noise parameter the smaller the steady state entanglement becomes. 
Parameters as in Figure 1. 

we begin from the "central" subspace (basis Wtg)^ \9^Y) of our density ma- 
trix, which for 7 becomes a decoherence free subspace Q. Wth these 
initial conditions we studied the Bell state ■i/'in = ^^^^^'^^ with matrix repre- 
sentation 
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(23) 



the initial mixed state 

ft„ = l/2|»)tol+l/2(^^^^iMKM±W)) (24) 

with matrix representation 
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(25) 



and the product state \eg) with matrix representation 



7 
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Figure 3: Linear entropy as a function of time and noise parameter M for Bell state 
ipi^ = l££>^££>. Parameters as in Figure 1. 
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We present our results for zero temperature in Figures 1 to 6. Figure 
1 shows the expected monotonous degradation of entanglement due to ex- 
ternal noise for Bell state ipin = ^^'^^^^'^ and product state ^/'m = \gg)- In 
Figure 2 we have the concurrence collapse and revival as a function of time 
and the noise parameter M for the Bell state i\)in = l^^l^^. We observe the 
expected collapse and revival of entanglement with the destructive effect of 
noise. But it is interesting that noise does not affect the concurrence mono- 
tonically as can also be seen in Figure 3 for the evolution of linear entropy. 
In Figure 4 we plot the collapse time and the mixedness as a function of 
the noise strength, to show this lack of monotonicity more explicitly. We 
call this negative effect '"stochastic antiresonance" '. In this Figure we also 
present the dependence of Manures as a function of the parameters A and 
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Figure 4: Sudden death time (upper left) and steady mixedness (upper right) as a function 
of noise parameter M for the Bell state Vin = Parameters as in Figure 1. The 

non monotonous dependence on M is obvious. Mantiresonance as a function of anisotropy 
A (lower left) and population relaxation rate 7 (lower right). We see that both A and 7 

have a small influence on Mantiresonance- 
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Figure 5: Concurrence collapse (left) and revival (right) as a function of time and noise 
parameter M for mixed state pi„ = 5 1 55) (55 1 + ^ ( '■l^^^"'"!^'^^-!^^'^^!"'"^^'^!-' ). Parameters as in 
Figure 1. 




Figure 6: Concurrence production and collapse (left) and concurrence revival (right) as 
a function of time and noise parameter M for the product state tpin = \eg)- We observe 
the stochastic antiresonance behavior in the entanglement oscillations of the left diagram. 
Parameters as in Figure 1. 
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7- The diagrams show that this dependences is non-zero but very smalL In 
Figure 5 we have the concurrence collapse and revival for the mixed state 
Pin = ^\99){99\ + l( '^l'^^^^l^'^^-'^^^'^^l^^^'^l^ ). In Figure 6 we have the concur- 
rence production and collapse and concurrence revival for the product state 
\eg). It is obvious that noise affects the entanglement production in an un- 
expected way. The entanglement oscillation amplitudes exhibit stochastic 
antiresonance behavior. Furthermore we see that for bigger values of M the 
oscillations last longer. 

5. Concurrence and purity evolution for T 7^ 

For T > the master equation ([T]) becomes: 

^ = - z[H, p,] + 7(n + l)D[S^]p, + 7(n + l)D[S^]p, 

+ jnD[St]ps + inD[S+]ps - M[V, [V,ps]], (27) 

where n denotes the average excitation quanta of the bath. It depends mono- 
tonically on the temperature and it is used to parametrize it. For T = we 
have n = and for T — )■ 00, n — > 00. The steady state concurrence is found 
to be equal to: 




(^73 + 4 fc2M72 + kn^-f + 2A^Mf 

' (fc2 {n + 1) n73 + 4M {{n^ + n) 72 + iA2) k + {Q^n{n + 1) + ^A^) 7)^ 
F (fcV + 4A;2M72 + ^27 + 2A2M)^ 

where k = {n + 1/2) and Q"^ = u"^ + A"^ + AM"^ 

For given system parameters Csteady decays rapidly when noise and ther- 
mal excitation parameters increase. Moreover n is a multiplicative factor in 
the Master equation and degrades the quantum entanglement. 

In Figure 7 we plot the concurrence collapse as a function of time and 
noise parameter M for three different values of the average photon number. 
We observe that the increase of temperature decreases the nonmonotonicity 
of the noise dependence, something normally expected. 

6. Conclusion 

We studied the entanglement evolution of a two-qubit Heisenberg XY 
chain in the presence of a noisy magnetic field. We solved the appropriate 
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Figure 7: Concurrence collapse as a function of time and noise parameter M for the Bell 
state ipin — l^^^^^''^ (J A = 0.1,7 — 0.01 and lo — 1), for average photon numbers 
n=0.02, n=0.06 and n=0.1, from left to right. The non monotonicity of the noise influence 
diminishes progressively. 
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Master equation both for T = and T 7^ and showed that our system 
presents stochastic antiresonance behavior for certain initial states. These 
states have matrix representation of the form of equation (20). E| 

For 7 — )• the system acquires a decoherence free subspace (DFS), 
spanned by the basis {\eg), \ge)}. Stochastic antiresonance in our system 
has to do only with the sudden death part of entanglement and not with the 
rebirth part, something that has advantages and disadvantages. 

Noise can play a destructive role in the first moments of evolution. There 
is a value of noise for which the system has the worst response (the small- 
est entanglement collapse time) and must be avoided, unlike the stochastic 
resonance case where we try to reach this value. This can be achieved easily 
because our noise source is classical and we assume that it is exactly tunable. 
The above have been confirmed both with the value of the concurrence and 
that of linear entropy. The steady state concurrence is affected negatively by 
the external noise. As for the temperature, it affects drastically the system 
evolution, diminishing the sharpness of antiresonance. 

The non monotonic dependence of the disentanglement time, of course 
makes sense only in the case of a finite time, that is in the case of Entangle- 



ment Sudden death (ESD) j23|, M, \M m \M |28|, [351 |36|, l3 



It is interesting 

m 



38, 39 



to investigate the influence of the geometry of the state manifold 
on this interplay of noise, environment and qubit interactions. We are cur- 
rently looking at these issues. 



^For initial states of this form, the first moments of entanglement evolution (sud- 
den death for initially entangled states and entanglement production for initially product 
states) are quantified by C2 while the rebirth-steady state entanglement is quantified by 
Ci. 
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7. Appendix 

7.1. Exact solution of the master equation 

The differential system that represents our master equation is: 

piiit) = -27(n + l)pii(t) + (pi4(t) - Pii{t))Ai + 7n(p22(t) + PssW) 

Pi4(t) = (-(2n - 1)7 - 2iu - AM)p^i{t) + (pn(t) - pAi{t))M 

P22it) = -fin + l)pii(t) - (2n + l)7P22(t) + Ji{p2z(t) - p^2(t)) + -fnpuit) 

P23{t) = -(2n + l)7P23(t) + Ji{p22{t) - P33{t)) 
P32(t) = -(2n + l)7/032(t) - J^(p22(t) - PSSW) 

P33(t) = 7(« + l)Pll(^) - (2n + l)7P33(t) + ^^(P32(t) - P23(t)) + 7^4^ 

P4i(t) = (-(2n + 1)7 + 2iu - AM)p^i{t) - {pi,{t) - pu{t))Ai 
PAiit) = 7(n + l)(p22(t) + P33(t)) + A2(p4i(t) - PmW) - 27np44(t) 

We observe that noise enters only in pi4 and p4i. Indeed, it can be easily 
shown that for 7 0, the system acquires a decoherence free subspace 
(DFS), spanned by the basis {\eg), \ge)}. This means that if the initial state 
of the system resides in the central subspace of the system Hilbert subspace 
in a symmetric fashion (Bell states ^^^^^^), then it remains the same. This 
central subspace of the general density matrix plays a crucial role for the 
unexpected resonance like behavior of certain initial states. 

The introduction of noise in the master equation makes the general so- 
lution quite complex. In the steady state limit the solution of our system 
is: 

1 (47n + 7 + mk) A2 + 477.27 (A;272 + AMk-f + Q^) 
16 (y^s^s + 4MA;272 + ^72^7 + 2MA2) k 

(2i7n + AiM + 77 + 2uj) 7A 
"SfcV + 32A:2M72 + m'^k-f + I6MA2 
(7 + SkM) A2 + (n2 + n) (4^273 ^ 4^ (^2 ^ /^^m^)) 
16A;27 (^272 + 4A;M7 + Vfi) + 32A;A2M 




P22(^) 

P44 = 1 - Pll - P22 - P33 

where k = {n + 1/2) and ^^2 = cj2 + A2 + 4M2. In this case C = Ci, because 
P32(^) — and C cannot take negative values. 



Pii(t) = 

P22{t) = 

P2z{t) = 
P32(t) = 
P33(^) = 
Pil{t) = 
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